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Recent work has demonstrated the importance of post-transcri- 
ptional gene regulation in toxic responses. In the present study, 
we used two rat models to investigate mRNA translation in 
the liver following xenobiotic-induced toxicity. By combining 
polysome profiling with genomic methodologies, we were able 
to assess global changes in hepatic mRNA translation. Dio3 
(iodothyronine deiodinase type III) was identified as a gene that 
exhibited specific translational repression and had a functional 
role in a number of relevant canonical pathways. Western blot 
analysis indicated that this repression led to reduced D3 (the 
protein expressed by Dio3) levels, enhanced over time and with 
increased dose. Using Northern blotting techniques and qRT-PCR 
(quantitative reverse transcription-PCR), we confirmed further 
that there was no reduction in Dio3 mRNA, suggesting that 



translational repression of Dio3 is an important determinant 
of the reduced D3 protein expression following liver damage. 
Finally, we show that drug-induced hepatotoxicity appears to 
cause localized disruptions in thyroid hormone levels in the liver 
and plasma. We suggest that this leads to reduced translation 
of Dio3 mRNA, which results in decreased D3 production. It 
may therefore be possible that this is an important mechanism 
by which the liver can, upon early signs of damage, act rapidly 
to maintain its own energy equilibrium, thereby avoiding global 
disruption of the hypothalamic-pituitary-thyroid axis. 

Key words: InnoMed PredTox consortium, iodothyronine 
deiodinase type III, liver, polysome profiling, toxicogenomics, 
translation. 
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INTRODUCTION 

The control of gene expression downstream of transcription is 
of physiological and toxicological importance due to the speed at 
which this level of regulation can be used to generate new proteins 
[1]. Such responses are essential within the CNS (central nervous 
system), for example, where cells demonstrate rapid variation 
in their metabolic activity [2,3], and within the liver where a 
rapid response is necessary to combat the constantly changing 
chemical milieu. Translational regulation enables a faster, more 
flexible response to cellular stress because the need to transport 
newly synthesized mRNAs to make more protein is avoided [3]. 
This response is demonstrated in cases of temperature shock 
and DNA damage where the majority of mRNA translation is 
shut down, with just a few key mRNAs increasing their rates of 
translation to enable the cell to deal with the additional stress [4] . 
Analysis of the extent to which mRNAs are recruited to ribosomes 
gives a quantifiable measure of translation efficiency and is 
a well-established technique [5]. This process can be coupled 
with genomics methodologies to allow a global assessment of 
the translational activity of mRNAs following cellular stress. 
Furthermore, in combining these data with those obtained from 
transcriptional studies, a comprehensive overview of the gene- 
regulatory processes used by the cell to adapt to, and recover 
from, stress can be obtained. This approach has demonstrated 
a role for both transcriptional and translational regulation in 
many pathophysiological states such as drug resistance and cell- 



cycle control [6,7]. In addition, with miRNAs (microRNAs) and 
other non-coding RNA species already established as translational 
regulators [8], research focusing on those regulatory processes 
occurring at the level of mRNA translation has increased. In fact, 
in certain cases, this level of control is the primary determinant 
of gene expression [2,9,10]. 

To investigate whether such translational mRNA regulation 
was important following drug-induced hepatotoxicity, we used 
two in vivo models. One of these was developed by the 
European Union Innovative IVIedicines Initiative (InnofVIed) 
PredTox consortium. The consortium undertook a series of 
experiments to investigate liver injury in the rat following 
dosing with novel pharmaceutical reagents that had failed during 
development due to overt toxicity [11]. The other was the well- 
established rat thioacetamide model [12,13], which we have 
included to supplement and verify the findings from the PredTox 
model. Using polysome fractionation and microarray methods, 
we were able to globally analyse differential mRNA translation, 
and through subsequent pathway analysis with Ingenuity IPA 
software (Ingenuity® Systems), we explored canonical pathways 
of potential interest and relevance in cases of hepatotoxicity. We 
identified Dio3 (iodothyronine deiodinase type III) as one of the 
genes that exhibited specific translational control under conditions 
of hepatic toxicity. 

Dio3 encodes an enzyme (D3) that is vital for TH (thyroid 
hormone) regulation; the maintenance of TH levels throughout 
life is of fundamental importance. The two major THs are T4 
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Abbreviations used: ALT, alanine transaminase; AST, aspartate aminotransferase; Cy3, indooarbocyanine; Cy5, indodicarbocyanine; D1-D3, 
iodottiyronine deiodinase type l-lll (encoded by Dio1-Dio3 respectively); GAPDH, glyceraldefiyde-3-pfiosptiate defiydrogenase; IcRNA, long non-coding 
RNA; qf^T-PCR, quantitative reverse transcription-PCR; RXR, retinoid X receptor; T2, 3,3'-diiodottiyronine; T3, 3,5,3'-triiodottiyronine; rTa, reverse T3; T4, 
ttiyroxine; TBST, Tris-buffered saline witti 0.05% Tween 20; TH, tfiyroid tiormone; TR, TH receptor. 

^ Correspondence may be addressed to either of tfiese autfiors (email kmp15@leicester.ac.uk or tim.gant@phe.gov.uk). 

The full dataset for translational microarray analysis of liver samples treated with PredTox hepatotoxic compound FP014SC has been submitted to the 
Gene Expression Omnibus under accession number GSE38807. 
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(thyroxine), secreted by the thyroid gland, and its biologically 
active form, T3 (3,5,3'-triiodothyrorune). In combination, these 
two molecules play a critical role during development, cellular 
proliferation and metabolic homoeostasis [14—16], and even 
transient disruptions in their levels can cause a wide number of 
transcriptional gene changes [17]. Prolonged deficiency causes 
growth retardation, impaired cognitive function and, in severe 
cases, infantile haemangioma [18,19]. In addition to D3, there 
are two further deiodinases (Dl and D2), all of which are highly 
conserved throughout the vertebrate kingdom [20]. Of the two 
THs, T4 is the most abundant; it is transported through the 
cellular membrane, via membrane transporters, before diffusion 
to the nucleus where it is converted into the more active T3 
by Dl and D2 [21]. On a tissue-specific basis, if levels of 
T3 and T4 become too high, then D3 is recruited from the 
plasma membrane [22], and, together with Dl, it converts T3 
and T4 into the inactive metabolites 1T3 (reverse T3) and T2 (3,3'- 
diiodothyronine) respectively. Levels of D3 are T3 -dependent, 
requiring the presence of TR (TH receptor) a [23]; this feedback 
mechanism is vital for the control of both the hypothalamic- 
pituitary-thyroid axis [24,25] and TH levels throughout the body. 

Dio3 shows particularly high expression in fetal and placental 
tissues [26,27], where it functions to protect the fetus from high 
levels of maternal TH. The majority of adult tissues demonstrate 
only low expression of Dio3; however, studies show that the rat 
brain maintains high D3 activity throughout life [28]. These data 
supports a functionally protective role for the enzyme, in this 
instance in modulating thyroid levels within neurons. In addition, 
in some circumstances, there is re-expression of Dio3 in the adult, 
for example during proliferation and cell growth [29]. This has 
been demonstrated in a number of pathophysiological conditions 
including cancer, myocardial infarction and liver regeneration 
following partial hepatectomy [19,30,31]. 

Although much is known about the transcriptional regulation 
of Dio3, its translational, or post-transcriptional, regulation has 
not been extensively investigated. Levels of Dio3 are increased 
transcriptionally in the presence of THs [32], retinoic acid [33] 
and growth factors [34], and reduced by growth hormones 
[35] and in hypothyroidism. Genomic imprinting of Dio3, via 
differentially methylated regions [36], is also suggested to be a 
contributor to transcriptional regulation [32]. Although not yet 
identified in the case of Dio3, there are published data on the 
other deiodinases that suggest that post-transcriptional regulation 
is essential; for example, Dio2 demonstrates significant post- 
transcriptional regulation [37]. Within the brain, the translational 
regulation of Dio2 is critical for the preservation of TH levels to 
prevent hypothyroidism, and it has been shown that the increases 
in D2 activity are far greater than can be accounted for by 
transcriptional changes in Dio2 mRNA [38]. Within the rat and 
human adult liver, Dio2 is not expressed [39]; however, Dio3 is. 
We therefore investigated whether within the liver this important 
mode of regulation could be controlled by Dio3 rather than Dio2. 

In the present study, we show that Dio3, like Dio2, can 
show a change in expression at the post-transcriptional level. 
We hypothesize that this mechanism may play a key role in 
the synthesis of D3 protein, the control of TH levels and the 
maintenance of the hypothalamic-pituitary- thyroid axis in cases 
of drug-induced cellular stress to the liver. 

EXPERIMENTAL 
Animals: PredTox 

With a view to reducing the financial and time restraints caused 
by the extensive testing required for any new drug candidate, a 



consortium comprising 15 pharmaceutical companies, two SMEs 
(smeill and medium enterprises) and three universities was set 
up in 2005. This consortium, the InnoMed PredTox project, 
formed under Eiu-opean Union Framework Progranune 6, sought 
to develop genomic biomarkers for the early detection of drug- 
induced toxicity. Male Wistar rats received pharmacologically 
active doses of 14 compounds that had failed in development 
owing to hepato- and/or nephro-toxicity. The treatment regime 
involved dosing for up to 15 days with a vehicle dose, a low dose 
or a high dose of each compound. The five agents that caused the 
most pronounced hepatobiliary injury were selected for further 
investigation in the present study. Liver samples from rats treated 
with these five compounds were provided as donations from the 
PredTox consortium. The compounds and their properties are 
shown in Supplementary Table SI at http://www.biochemj.org/bj/ 
453/bj4530071add.htm. 

Animals: tiiioacetamide model 

Male Wistar rats (3 months, 300-350 g) were dosed via the 
i.p. (intraperitoneal) route with Thioacetamide (Sigma) dissolved 
in 0.9% saline to final doses of 50 mg/kg, lOOmg/kg and 
150 mg/kg. These doses fall below the threshold for overt toxicity 
[12,13]. Vehicle-only-treated rats were used as controls. Six 
animals were used at each dose level and all Figures show 
means + S.E.M. for each measurement. After 24 h, half of the 
animals from each group were anaesthetized using isoflurane and 
blood was withdrawn via the descending vena cava. Blood was 
collected into lithium/heparin tubes and the plasma was isolated. 
Organs were perfused in situ with PBS and animals were killed by 
decapitation under terminal anaesthesia. Organs were harvested, 
sections were taken for histological analysis and the remainder 
were snap-frozen in liquid nitrogen. The remaining animals from 
each group underwent the same process after a further 24 h. All 
procedures were licensed under U.K. Home Office project licence 
80/2126. 



IVIeasurement of ALT (alanine transaminase) and AST (aspartate 
aminotransferase) levels 

To measure plasma levels of ALT and AST, kits from Sentinel 
Diagnostics (Alpha Laboratories) were used according to the 
manufactiu-er's instructions. 



Translational profiling 

Sucrose density fractionation was used to separate mRNAs 
according to their ribosomal mass. mRNAs that are under active 
translation will have a large number of ribosomes attached 
[40], increasing their density and resulting in them migrating 
towards the bottom of the gradient; those that are less actively 
translated, with fewer ribosomes attached, will be less dense 
and will therefore be situated towards the top of the gradient. 
Approximately 150 mg of each liver sample was ground to a 
powder in liquid nitrogen using a pestle and mortar and then lysed 
in 15 mM Tris/HCl (pH 8.0), 300 mM NaCl and 15 mM MgClj, 
plus inhibitors [1 mg/ml heparin, 100 /xg/ml cycloheximide 
and 80 units of RNasin® (Promega)]. Following centrifugation 
at 12 000 g for 5min at 4°C, the resulting supernatant was 
added to a 10-50% sucrose gradient and ultracentrifuged at 
388000 rev./min in a Beckman ultracentrifuge for 2h using a 
SW40 rotor. Sucrose of a greater density was used in combination 
with a Harvard syringe pump system to divide each gradient into 
ten to 16 individual 1 ml fractions. These fractions were collected 
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entire mouse genome. We had established previously that these 
microarrays showed excellent cross-species reactivity between 
mouse and rat (T. Gant, unpubUshed work). Following overnight 
incubation, the sUdes were washed [wash 1: Ix SSC (Ix 
SSC is 0.15 M NaCl/0.015 M sodium citrate) and 0.03 % SDS 
for 5min; wash 2: 0.2x SSC for 3 min; and wash 3: 0.05 x 
SSC for 3 min] and scanned on a 4200A Axon scanner 
(Molecular Devices). Experiments were performed in duplicate, 
incorporating a dye-swap technique for n = 5 pairs of samples. 
Once all hybridization reactions were complete, the results 
were normaUzed by LOWESS (locally weighted scatterplot 
smoothing), which uses locally weighted regression to smooth 
scattered data (S. Zhang, personal communication) and tested 
for statistical significance using a two-tailed Student's t test as 
described previously [41]. There were a maximum of four values 
for each mRNA, corresponding to the proportional representation 
of it within each subpool of fractions. By calculating the change 
in values, i.e. degree of slope, across the monosomal, light 
polysomal, medium polysomal and dense polysomal region, it 
was possible to determine any translational shift. In addition, the 
overall change in transcriptional activity across the gradient was 
calculated for each mRNA. The full dataset has been submitted 
to GEO (Gene Expression Omnibus) under accession number 
GSE38807. Pathway analysis of differentially translated mRNAs 
was performed using Ingenuity IPA software. 



Figure 1 Translational profiling of iiver samples talcen from rats treated 
witli tlie PredTox compound FP014SC 

RNA from the livers of control (vehicle-treated) and test (compound-treated) animals was 
fractionated by sucrose density centrifugation and pooled into four fractions: monosome (M), 
light polysome (LP), medium polysome (MP), and heavy polysome (ffP). Each of the four control 
(C) fractions were labelled with Cy3 and hybridized against each of the four corresponding test 
(T) fractions labelled with Cy5. The statistical significance of differences between C and T in 
each of the four fractions was calculated using a reverse-labelled two-tailed Student's I test 
[41]. All data with four P values >0.05 were rejected. Linear regression analysis was then 
performed on the four log2 ratios of test fractions/contral fractions. This showed movement 
of mRNAs through the gradient and enabled differentially translated mRNAs to be identified. 
Reverse-labelling reactions were also performed and incorporated into the analysis. 

into tubes containing 1 ml of Tri Reagent (Sigma), and RNA was 
extracted according to the manufacturer's instructions. 



Genoinic analysis of mRNA fractions 

The individual fractions were subpooled into four groups, 
monosomal, light polysomal, medium polysomal and heavy 
polysomal (Figure 1), according to ribosomal density. Microarray 
analysis was performed by hybridizing control (vehicle-treated) 
monosomes against test (high-dose-treated) monosomes on one 
microarray, control light polysomes against test light polysomes 
on a second microarray, and so forth for each subpool of fractions 
(Figure 1). Equal amoimts of RNA from each subpool were 
precipitated with 10 /xl of 3 M sodium acetate (pH 5.2) and 275 /xl 
of 100% ethanol at — 20°C overnight. Following centrifugation 
at 12000^ for 10 min at 4°C, pellets were washed twice with 
75 % ethanol and resuspended in \Q [A of water. The precipitated 
RNA was reverse-transcribed, hydrolysed and coupled to a dUTP- 
conjugated Cy3 (indocarbocyanine) or Cy5 (indodicarbocyanine) 
dye (GE Healthcare). The labelled cDNA samples were mixed 
and hybridized at 42 °C overnight to a 70-mer MEEBO (Mouse 
Exonic Evidence-Based Oligonucleotide) microarray (http:// 
www.microarray.org/sfgf/meebo.do). The microarrays were 
printed in-house using an Array Jet Ultra Marathon microarrayer 
and an lUumina probe set (Invitrogen), which covers the 



Western blotting 

Tissue samples were lysed at 4°C in lysis buffer [50 mM 
Tris/HCl (pH 7.4), 150 mM NaCl, 5 mM EDTA, I % Nonidet 
P40, 0.25% sodium deoxycholate and protease inhibitors]. 
Protein (30 /ig) was separated on a 10% acrylamide gel and 
transferred on to a nitrocellulose membrane (GE Healthcare). 
Membranes were blocked with 10% (w/v) Marvel non-fat dried 
skimmed milk powder in TBST (Tris-buffered saline with 0.05 % 
Tween 20) and incubated at 4°C overnight with the primary 
antibody anti-Dio3 (1:1000 dilution) (Novus Biologicals) in 5 % 
Marvel. This antibody has previously been validated [42,43] 
against the only other available anti-D3 antibody [44]. Following 
overnight incubation and washing in TBST, membranes 
were incubated with the appropriate horseradish peroxidase- 
conjugated anti-rabbit secondary antibody (1:1000 dilution), 
(Santa Cruz Biotechnology) and visualized using ECL (enhanced 
chemiluminescence) Western blotting detection reagents (GE 
Healthcare). To check for equal loading, membranes were stripped 
using 62.5 mM Tris/HCl (pH 6.8), 2 % (w/v) SDS and 100 mM 
2-mercaptoethanol, re-blocked and incubated overnight with anti- 
GAPDH (glyceraldehyde-3-phosphate dehydrogenase) (1:30000 
dilution) (Sigma). Western blots were quantified using the gel 
analysis function of Image J software (NIH). To analyse time- 
course and dose-response effects, one-way ANOVA with Tukey's 
multiple comparison post-hoc test was performed using GraphPad 
Prism (version 5.01). Results are expressed as means + S.E.M. 
P < 0.05 was considered significant. 



RNA isolation and qRT-PCR (quantitative reverse 
transcription-PCR) 

RNA was extracted from <100mg of Uver samples using 
Tri Reagent according to the manufacturer's guidelines. RNA 

quality was assessed using an Agilent Bioanalyser, where a 
RIN (RNA integrity number) score of >8.5 was used as an 
indication that RNA was of a suitable quality for further work. 
To prevent genomic DNA contamination, 10 /xg of RNA was 
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Table 1 Summary of the liistopatliological data for eacli of tlie PredTox compounds 

Severity ievels: +++>++>+. 



Frarreworl< Programme 6 study name Compound sponsor 



Toxicological response within iiver 



Severity 



FP004BA 



FP005IVIE 



FP007SE 



FP013N0 
FP014SC 



Bayer 



Merck 



Boeliringer Ingelheim 



Novartis 



Bayer Schering Ptiarma (tormerly Schering) 



Bile duct damage 
Hepatocyte necrosis 
Regenerative hyperplasia 
Intlammatory responses 
Cholestasis 
Fibrosis 

Hepatocellular apoptosis 
Necrosis 

Peribiliary intlammation 
Fibrosis 

Bile duct proliferation and necrosis 
Hepatocellular hypertrophy 
Increased transamininases, ALP and bilirubin 
Pericholangitis with bile duct hyperplasia 

Cholestasis 
Intlammation 

Hepatocellular hypertrophy 
Vacuolation 

Absence of hepatocellular glycogen deposits 
Increased tafty deposits 
Severe acute necrotic liver injury 
Elevated liver transaminases 
Cholestasis 

Hepatocellular vacuolation and hypertrophy 
Regeneration and Increased mitosis 



++ 



+++ 



treated with Turbo DNase (Applied Biosystems) according to the 
manufacturer's instructions. The DNase was inactivated tlirough 
re-extraction of the RNA, as described above. RNA (1 /ig) was 
reverse-transcribed using a high-capacity reverse transcription kit 
(Applied Biosystems) and then qRT-PCR was performed using an 
ABI Prism® 7000 Sequence Detection System. Primers against 
Dio3 and Actb (/S-actin) (housekeeping gene) were designed 
using Primer Express version 2.0.0 (Applied Biosystems) and 
ordered from Sigma-Aldrich. Primer sequences are given in 
Supplementary Table S2 at http://www.biochemj.org/bj/453/ 
bj4530071add.htm. All primers were tested for percentage effi- 
ciency before use [45]. Reactions were performed in 25 /xl of 
Power SYBR® Green PGR Master Mix (Applied Biosystems), 
using 450 nM each of forward and reverse primer. 

Northern blot analysis 

An equal volume of RNA from each fraction was loaded 
on to a denaturing agarose gel as described previously [46] 
and, following overnight electrophoresis, was transferred on to 
MagnaGraph nylon transfer membrane (Labtech International) 
by capillary transfer (lOx SSC). The membrane was cross-linked 
with 120 mJ/cm^ UV for 60 s. In addition, fractionated RNA was 
mixed with denaturing buffer [59 % (v/v) deionized formamide, 
24% (v/v) formaldehyde and 0.14 M Mops], dot-blotted on to 
the same membrane and cross-linked as described above. Probes 
derived from PGR products to the coding regions of Dio3 and 
Actb (primer sequences given in Supplementary Table S2) were 
labelled with '^P using Escherichia coli Klenow fragment, and 
incorporation was measured by scintillation counting. Probe (10* 
d.p.m./ml) was added to the pre-hybridized [50 % (v/v) deionized 
formamide, 6x SSG, 5x Denhardt's solution (0.1% Ficoll 
400/0. 1 % polyvinylpyrrolidone/0. 1 % BSA) and 1 % (w/v) SDS] 
membranes and hybridized for 48 h at 42 °G. The hybridized 
membranes were visualized using autoradiography and quantified 



using ImageQuant TL (version 2003.03). The value for each 
gradient fraction was expressed as a percentage of the total for 
the whole gradient. 



Determination of T3 and T4 concentrations 

The levels of T3 and T4 in the liver and plasma were determined 
by highly sensitive and specific RIAs, as described in detail in 
[47,48]. 



RESULTS 

Assessment of liver damage 

Previous work by the PredTox consortium indicated that the 
five compounds selected for further analysis in the present 
study caused hepatic injury, as determined by histopatiiology 
and clinical chemistry (Table 1). We supplemented the PredTox 
model with the thioacetamide-treated rat model. To determine the 
degree of liver injury in this additional model, plasma ALT and 
AST levels were measured 24 h after administration of 100 mg/kg 
thioacetamide. Enzyme levels were increased more than 10 (ALT) 
and 43 (AST) -fold when compared with the vehicle-treated 
samples (Figure 2A). Levels of both enzymes in the treated 
animals were lower at 48 h than at 24 h, reflecting the acute 
nature of the hepatic injury and subsequent initiation of repair 
mechanisms. 

Glinical chemistry data were supported by histopathological 
analysis. Representative images from the livers of rats treated with 
the vehicle and the highest dose (150 mg/kg) of thioacetamide 
are shown in Figure 2(B). The images from the treated livers 
demonstrated some fibrosis and inflammation; the hepatocytes 
had small dense nuclei, indicative of hepatic injury, and there was 
also evidence of hepatocyte vacuolation (Figure 2B, indicated 
by V). However, there were also signs of some regeneration 
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Figure 2 Clinical chemistry and histopatliology indicate liepatotoxicity foilQwing treatment witli tliioacetamide and tiie PredTox compounds 

(A) Plasma ALT and AST levels (in units/litre) and (B) histopathologlcal images 48 h after treatment with (1) vetiicle-only and (11) 150 mg/kg ttiioacetamide. V, vacuolation, indicative of fiepatocellular 
injury; R, area rich in hepatocytes forming parenchymal nodules, indicative of enhanced regenerative activity. (C) ALT and AST levels relative to control animals from the plasma of rats treated with 
the PredTox compound FP014SC. Results are means + S.E.M. (n = 3). *P ^ 0.05, "P ^ 0.01 , *"P < 0.001 (ANOVA with Dunnett's post-hoc test to compare groups). 



(Figure 2B, indicated by R) after 48 h, concurrent with the 
decreased levels of serum transaminases. 

ALT and AST levels were similarly increased following 
treatment with the PredTox compounds. Representative data 
from the most severe hepatotoxic agent (FP014SC) is given in 
Figure 2(C). As with the thioacetamide, this effect was transient 
and levels had returned to close to those seen in the vehicle-treated 
samples by 15 days after treatment. 

The histopathologlcal summary for each of the PredTox 
compounds indicated a histology similar to that seen in the 
livers of the thioacetamide-treated rats, with common features 
including apoptosis, inflammation, cholestasis and regenerative 
repair (Table 1). 

MIcroarray analysis indicated that D3 Is Involved In drug-induced 
hepatotoxiclty 

Tremslational microarray analysis was performed on RNA from 
those PredTox samples treated with the most potent hepatotoxic 
agent, FP014SC (Table 1), at the highest dose and latest time point 
(1120 mg/kg, day 15). Analysis of polysome-associated mRNAs 
showed an alteration in the numbers of ribosomes recruited to indi- 
vidual mRNAs, indicative of changes in translational efficiency. 
All mRNAs that demonstrated a shift in ribosomal occupancy 
across the gradient (Supplementary Table S3 at http://www. 
biochemj.org/bj/453^j4530071add.htm) were subjected to path- 
way analysis using IPA software. From this analysis, 169 
canonical pathways were identified as having two or more altered 
genes (Supplementary Table S4 at http://www.biochemj.org/bj/ 
453/bj4530071add.htm). Pathways under the control of RXR 
(retinoid X receptor) were highlighted because RXR plays 



a fundamental role in many disrupted metabolic pathways. 
Individual examination of RXR-based canonical pathways 
identified Dio3 as an indirect regulator of RXR, via the control 
it exerts over TH levels. This is demonstrated in the TR/RXR 
activation pathway (Figure 3); other related canonical pathways 
(indicated by CP) are also highlighted. Dio3 displayed a shift 
from the heavy polysomal region, to the lighter monosomal 
region following treatment with FP014SC (Supplementary Table 
S3), reflecting decreased ribosomal occupancy. This indicated 
that Dio3 was less actively translated following treatment with 
FP014SC. 



D3 Is down-regulated following drug-Induced hepatotoxiclty 

Reduced translational efficiency of an mRNA should cause 
a corresponding decrease in protein levels. We carried out 
Western blot analysis to determine D3 levels in the livers of 
rats treated with the PredTox compounds or thioacetamide. 
Following high-dose treatment for 15 days with FP004BA, 
FP007SE and FP014SC (Figures 4A-4C), there was a significant 
reduction in D3. In addition, high-dose treatment (150 mg/kg) 
with thioacetamide for 48 h (Figure 4D) caused significantly 
decreased D3 levels. A similar trend was seen with the other 
PredTox compounds, FP005ME and FP013NO (Figures 4E 
and 4F), although they failed to reach statistical significance. 
When considering samples on an individual basis, there was 
a correlation between the extent of liver damage, according to 
clinical and histopathologlcal grading, and reduction in D3 levels. 
We therefore also tested whether dose- and time-course-related 
responses were evident. To measure time-related effects, we used 
the PredTox sample that had shown the largest reduction in 
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Figure 3 Networks of tlie TR/RXR activation patliway 

All mRNAs that demonstrated a ctiange In translatlonal eftlclency were subject to pattiway analysis using Ingenuity® IPA software. TR/RXR was Identified as a key pattiway. Associated RXR canonical 
pathways are also highlighted ('CP'): FXR, farnesold X receptor; IL-1, Interleukin 1; LPS, llpopolysaccharlde; LXR, liver X receptor; PPARa, peroxisome-proliferator-activated receptor a\ PXR, 
pregnane X receptor Dio3 plays a regulatory role In the TR/RXR pathway via the control It exerts over thyroid levels. 



D3 levels (FP007SE) following high-dose treatment for 15 days. 
The reduction in D3 was exacerbated over time, until, by day 
15, levels were reduced to less than 15% of the vehicle- 
treated samples (Figure 5A). We used the thioacetamide model 
to look for a dose-response effect in the reduction of D3 
with hepatotoxicity. Although the only statistically significant 
reduction in D3 was seen following treatment at the highest 
dose (150mg/kg), there was a small, but consistent, change 
at the lower doses (Figure 5B). Two of the other PredTox 
compounds caused a similar effect (Supplementary Figure SI 
at http://www.biochemj .org/bj/453/bj453007 1 add.htm). 



The reduction in protein was not due to ciianges in tlie levei of Dio3 
mRNA 

To confirm that changes in D3 at the protein level were 
occurring independently of a change in transcriptional activity, 
we performed qRT-PCR on total RNA. This showed either no 
change in Dio3 mRNA levels or, in the case of FP007SE, a 
statistically significant increase in Dio3 mRNA levels (Figure 6A). 
High-dose thioacetamide treatment caused a non-significant 
increase 48 h after dosing (Figure 6B) and the other PredTox 
compounds resulted in no change in Dio3 mRNA levels 
(Supplementary Figure S2 at http://www.biochemj.org/bj/453/ 
bj4530071add.htm). 



Less D3 was recruited to tlie polysomes following drug-induced 
liepatotoxicity 

Using the FP007SE and thioacetamide studies, we performed 
Northern blot analysis to verify that changes were occurring post- 



transcriptionally, at the translatlonal level. Hybridizing with a 
radiolabelled probe enabled the proportional representation of 
Dw3 across the density gradient to be determined. We saw an 
increase in Dio3 mRNA abundance in the lighter fractions of 

the compound-treated samples compared with the vehicle-treated 
samples (Figures 6C, i, and 6D, i), indicating lower translatlonal 
efficiency. A similar validation technique was performed by 
machine dot-blotting RNA from each fraction across the gradient 
on to a nylon membrane and hybridization with the same 
D;oJ-specific probe. As with the initial Northern blot analysis, 
when quantified, the treated samples demonstrated a shift from 
the heavier density fractions to the lighter density fractions 
(Figure 6C, ii, and 6D, ii). This supported the hypothesis that 
the changes in protein levels were due to fewer ribosomes being 
recruited and therefore less efficient translation of Dio3 mRNA. 



The liver and plasma had reduced levels of T3 and T4 following 
toxic injury 

We measured levels of T3 and T4 within the liver following 
treatment with FP007SE or thioacetamide. We saw that, following 
treatment with FP007SE for 15 days, there was a reduction in 
levels of both hormones: by 57 % (T3) and 53 % (T4) (Figure 7A). 
Similarly, thioacetamide caused a reduction in T3 and T4 levels 
by 73 % and 25 % respectively 24 h after dosing (Figure 7B, i). 
The levels of T4 were also reduced in two of the other PredTox 
studies (Supplementary Figure S3 at http://www.biochemj.org/bj/ 
453/bj4530071add.htm). We also measured the levels of T3 and 
T4 in the plasma of the thioacetamide-treated rats. T3 levels were 
reduced by up to 60 % 24 h after dosing. Although statistical 
significance was not achieved, T4 levels also showed a trend for 
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Figure 4 D3 protein is down-regulated following drug-Induced liver damage 

Western blot analysis was performed on lysates from the livers of rats treated at tiigh-dose levels with various fiepatotoxic agents. D3 protein levels were normalized to ttiose of GAPDH. Representative 
gel images are shown for PredTox compounds FP004BA (A), FP007SE (B), FP014SC (C) and thioacetamide (D). Non-significant reductions were seen v^ith PredTox compounds FP005ME (E) and 
FP013NO (F). IVlean relative levels of D3 are indicated beneath each gel image. The position of a 37 kDa protein is indicated. Results are means ± S.E.M. (n^3). *P ^ 0.05, **P ^ 0.001 (Student's 
ftest). 



a reduction following treatment (Figure 7B, ii). T4 has a long 
half-life in the plasma and therefore the concentration of this 
enzyme within the systemic circulation usually remains relatively 



stable [20]. Sample availability from the PredTox studies was 
limited, and did not allow measurement of T3 and T4 levels in the 
plasma. 
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Figure 5 D3 protein is down-regulated in botli a time- and dose-dependent 
manner 

Western blot analysis was performed on lysates from tfie livers of rats treated with various 
hepatotoxic agents. D3 protein levels v^ere normalized to those of GAPDH. Representative gel 
images are shown for (A) PredTox compound FP007SE and (B) thioacetamide. Mean relative 
levels of D3 are indicated beneath each gel image. The position of a 37 kDa protein is indicated. 
Results are means + S.E.M. (n=3). *P^0.05, **P^0.01, 0.001 (ANOVA with 
Dunnett's post-hoc test to compare groups). 



DISCUSSION 

In the present study, we used global genomic methods to 
investigate changes in the rate of translation of thousands of 
individual mRNAs concurrently in two models of drug-induced 
hepatotoxicity. A change in the rate of mRNA translation can 
bring about a rapid efficient change in protein levels, without the 
need to generate more transcript [3]. We identified Dio3 as an 
mRNA showing a reduction in the rate of translational activity 
following hepatotoxicity. The present study provides evidence 
that the translational repression of Dio3 results in less D3 protein, 
independent of any change in mRNA level, and suggests that 
the reduced translation is a regulatory mechanism, enabling the 
liver to combat cellular stress. The known function of Dio3 
in modulating thyroxine levels supports further the hypothesis 



that post-transcriptional regulation may be a rapid compensatory 
response to drug-induced damage. 

In the rat, D3 activity is reduced in cases of hypothyroidism 
[28,49], although the extent of the change varies between tissues 
[20]. The observed reductions in T3 and T4, along with the tran- 
scriptional down-regulation of Diol mRNA (Supplementary Fig- 
ure S4 at http://www.biochemj.org/bj/453/bj4530071add.htm), 
which is known to be down-regulated in cases of hypothyroidism, 
indicated that the livers of the treated rats were in a hypothyroid 
state [17,50]. 

Previously, changes in D3 activity have been found to correlate 
strongly with changes in expression of Dio3 mRNA [49] ; however, 
we found no evidence for the change in protein production 
being regulated at a transcriptional level. Moreover, one of the 
PredTox compounds (FP007SE) showed a transcriptional increase 
in Dio3 mRNA within the liver following high-dose treatment 
(Figure 6A). This was in agreement with previous studies on 
Dio3, where an increase in expression was investigated [29-31]. 
To our knowledge, all previous work on Dio3 in the liver has 
focused on cases of enhanced cellular proliferation, either during 
growth and development in the fetus [26] or in the adult at times 
of cellular stress, such as during critical illness [29]. Previous 
studies have assumed that the adult liver contains low to negligible 
levels of Dio3 mRNA [23,50] and that, to measure transcriptional 
changes, cellular proliferation must first be induced. The PredTox 
consortium reported that, following 15 days of treatment with the 
FP007SE compound, there were signs of liver recovery and regen- 
eration [51]. Furthermore, 48 h after dosing with thioacetamide, 
there were chnical and histological signs of liver repair, which 
concurred with an increase in Dio3 mRNA levels. Therefore the 
results of the present study support what is widely reported in the 
literature whereby, in cases of cellular proliferation, Dio3 mRNA 
shows an increase in expression. However, despite this increase 
we saw an almost total abolition of conversion of the mRNA into 
protein. We therefore hypothesized that the changes in protein 
levels were mediated predominantly by altered Dio3 translation. 

It is well established that another of the deiodinase genes, 
Dio2, is regulated at levels distinct from transcription and 
that these alter the levels of D2 expression independently of 
changes in mRNA level [37,52]. The change in activity of D2 
is often greater than the change in Dio2 mRNA, for example in 
brown adipose tissue following cold exposure [53] and the brain 
following experimentally induced hypo- and hyper-thyroidism 
[54]. This is essential for homoeostasis of the THs. Certainly, 
in cases of endoplasmic reticulum stress, D2 activity is reduced 
independendy of transcriptional changes and this leads to a rapid, 
but significant, decrease in the levels of T3 [37]. 

Taken in combination, the post-transcriptional control of D2 
is therefore of critical importance in maintaining TH levels in 
mammals. However, within the rat and human liver, D2 is not 
expressed [20], and, as a result, an alternative method of regulation 
is necessary. We hypothesize that, following hepatotoxicity, 
D3 is recruited for TH homoeostasis in the liver. This seleno- 
enzyme functions primarily to inactivate T3 and T4 by conversion 
into their inactive metabolites, rT, and T2 respectively [20]. As an 
inactivating enzyme, D3 acts primarily to protect tissues from an 
excess of TH. If this mechanism is disrupted in some way, severe 
hypothyroidism [19] or hyperthyroidism can occur. 

We propose a mechanism whereby the hepatotoxins used in 
the present study have led to a need for increased energy for 
cellular proliferation, as part of the inflammatory response within 
the liver following cell damage. This is achieved through a rapid 
decrease in Dio3 translation and consequent reduction in the level 
of D3 protein. Although we cannot rule out a role for increased 
protein degradation, the reduced ribosomal association of Dio3 
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Figure 6 The reduction in D3 protein is due to changes in mRNA translation not transcription 

qRT-PCR analysis was performed on RNA extracted from ttie livers of rats treated witfi high doses of (A) PredTox compound FP007SE for 15 days or (B) thioacetamide. Results are means + S.E.M. 
(n = 3). *P < 0.001 (Student's f test). Polysomal association of Dio3 was measured by Northern blot analysis. Liver samples from rats treated with high doses of either (C) FP007SE (1 5 days) or 
(D) thioacetamide (48 h) were lysed and fractionated using a 10-50% sucrose gradient. RNA was extracted from individual fractions and eiecfrophoresed on a denaturing agarose gel, before being 
transferred on to a nylon membrane (1). Sizes are indicated in kb. RNA from each fraction was also dot-blotted on to a nylon membrane (11). All membranes were probed with a ^^P labelled probe 
specific to Dio3. The proportional representation of Dio3 within each fraction was calculated. Representative images are sfiown for both compounds and results are ttie mean + S.E.M. (n=3) 
amount of DioS within each fraction as a percentage of the total (iii). 



mRNA shown experimentally in the present study in both model 
systems provides a strong case for translational repression being 
an important mechanism leading to the reduction in D3 protein. 
In cases where DioS transcription is significantly up-regulated, as 
seen in Figure 6, in addition to a reduction in mRNA translation, 
protein degradation could also contribute to the reduced 
D3 protein levels. 

The benefit of targeting Dio3 post-transcriptionally is that the 
required response is efficient and rapid. A reduction in protein 
levels means less of the inactive metabolites rTj and T2 are 
produced and alterations in T3 and T4 levels can be restored to 
normal (Figure 8). This is vital for maintaining the hypothalamic- 
pituitary-thyroid axis and means that all control is carried out 
at the local level, in concordance with the current literature [55]. 
As the liver regenerates and clinical signs of damage are reversed, 
protein levels are reduced further. This is probably for the required 
homoeostasis as the liver expends sufficient energy into tissue 
repair. 



Although the levels of T3 and T4 within tissues are largely 
independent of those seen in the plasma [56], and the expression 
of the deiodinases is generally under localized regulation, the 
activity rate of the three deiodinases does directly affect the levels 
of circulating T3 and T4. This enables the feedback loop between 
the systemic organs and the thyroid gland to be maintained and 
has been experimentally vahdated in the case of D2, as described 
previously [37]. In the present study, we have demonstrated 
that, following treatment with thioacetamide, a known inducer 
of hypothyroidism, the levels of T3 and T4 were reduced in the 
plasma. The pattern of this change matched the changes in THs 
seen in the Uvers of treated animals (Figure 7B), although, as 
predicted by the Uterature [20], the change in T4 was smaller. 
D3 is located within the plasma membrane, although the bulk of 
the molecule is extracellular [22]. This gives it ready access to, 
and control over, the levels of circulating THs. We propose that, 
because the liver is in rapid equilibrium with the plasma, it detects 
a change in THs and acts rapidly to prevent further production of 
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Figure 7 T3 and T4 levels within the liver and plasma are reduced following 
drug-induced liver damage 

Liver sarrples were taken for measurement using highly sensitive and specific RIAs [47,48] 
following high-dose treatment with (A) PredTox compound FP007SE and (B, i) thioacetamide. (B, 
11) Plasma samples were taken for measurement as above following treatment with thioacetamide. 
Results are means + S.E.M. (n > 3). *P ^ 0.05, **P ^ 0.01 , ***P ^ 0.001 (Student's ( test). 



inactive metabolites (Figure 8). Through negative feedback, this 
mechanism enables levels of T3 and T4 in the liver to return to 
normal and ultimately restores homoeostasis within the liver. 

Despite it being over 15 years since Dio3 was first cloned 
[57], the regulation of this gene is complex and not yet fully 
understood. The gene is imprinted and preferentially expressed 
from the paternal allele [58]. In the mouse and humans, a non- 
coding gene transcribed in the antisense orientation, DioSos, 
has been identified [59]; and very recent work indicates that a 
homologue of this gene may also be present in the rat [60] . DioSos 
is proposed to belong to the family of IncRNAs (long non-coding 
RNAs), which, although commonly found in the mammalian 
genome, are poorly understood [61]. It has long been known 
that the stability of TRs can be moderated by naturally occurring 
antisense RNAs [62] and now that IncRNAs have been shown to 
directly regulate ribosomal association of their coding equivalent 
mRNAs [63], there is the distinct possibility that Dio3 is post- 
transcriptionally regulated in this manner. In addition to this in the 
better characterized mouse and human orthologues, DioSos and 
DioS are partially overlapping and demonstrate inverse correlation 
of expression [59]. 

We conclude that the reduction in D3 expression levels during 
hepatotoxicity can be correlated with a change in DioS mRNA 
translation rates. We hypothesize that this is an example of post- 
transcriptional control, which may be regulated by an IncRNA, 
DioSos. We demonstrate that the magnitude of the decrease in 
protein is directly proportional to the extent of liver damage. 
This is probably due to the energy demand required for active 
repair processes to be initiated. TH homoeostasis has already been 
shown to play a role in a range of pathophysiological conditions 
and the results of the present study suggest that the pathway is also 
perturbed by drug-induced hepatotoxicity. We are able to provide 
evidence that this forms part of a rapid response by the liver and 
speculate that this occurs so that the initiation of a localized event 
remains consigned to the liver and does not affect the whole of 
the hypothalamic-pituitary-thyroid axis. 
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Figure 8 Proposed mechanism of action following drug-induced hepatotoxicity 

Steps in bold type are those experimentally demonstrated in the course of the present study Non-bold type is used for the hypothesized mechanism of action and closure of the feedback loop. 
During hepatotoxicity the liver is in a hypothyroid state and both the plasma and liver have reduced levels of T3 and T4. Consequently the expression of D3 protein is reduced so that fewer inactive 
metabolites of T3 and T4 are produced and levels of both can return to normal, both locally (liver) and via plasma signals systemically 
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Figure S1 D3 protein is down-regulated in both a time- and dose-dependent 
manner 

Western blot analysis was performed on lysates from the livers ot rats treated with various 
hepatotoxic agents. D3 protein levels were normalized to those ot GAPDH. Representative gel 
images are shown tor PredTox compounds FP004BA (A) and FP01 4SC (B). Mean relative levels 
of D3 are indicated beneath each gel image. Results are means + S.E.M. (n = 3). *P < 0.05, 
**P < 0.01 , ***P < 0.001 (ANOVA with Tukey's/Dunnett's post-hoc test). 



^ Correspondence may be addressed to either of ttiese authors (email kmp15@leicester.ac.uk or tim.gant@phe.gov.uk). 

The full dataset for translational mioroarray analysis of liver samples treated with PredTox hepatotoxic compound FP014SC has been submitted to the 
Gene Expression Omnibus under accession number GSE38807. 
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Figure S2 The reduction in D3 protein is not due to clianges in mRNA transcription 

qRT-PCR analysis was performed on RNA extracted from tlie livers of rats treated witti high doses of PredTox compounds FP004BA (A), FP005ME (B), FP013MO (C) and FP014SC (D) for 15 days. 
Results are means ± S.E.M. (n - 3). 




I 

Inl 



dose 




High 



Figure S3 T4 levels within the liver are significantly reduced following 
drug-induced liver damage 

Liver samples were taken for measurement using highly sensitive and specific RIAs following 
high-dose treatment with PredTox compounds FP004BA (A) and FP014SC (B). Results represent 
means ± S.E.M. (n ^ 3). *P ^ 0.05, **P ^ 0.01 (Student's f test). 
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Figure S4 Dio1 mRNA is transcriptionally down-regulated 

qRT-PCR analysis was performed on RNA extracted from the livers of rats treated with higti doses of PredTox compounds FP004BA(A) and FP007SE (B) for 15 days and (C) 150 mg/kg thioacetamide. 
Results are means + S.E.IVI. (n = 3). *P ^ 0.05, "P ^ 0.0C1 (Student's f test). 
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Table SI PredTox compounds and their properties 

All available data on the five PredTox compounds selected for furttier analysis in ttie present study are provided. 



Framework 

Programme 6 Compound 

study name sponsor Official name 



Proposed iiigti Proposed low 
Structure dose{mg/kg) dose(mg/kg) Target organs Ctiemical class Ptiarmacological class 



FP004BA Bayer ( + )-(1 fi)-1 -[4-(4-fluorophenyl)-2,6-di-isopropyl-5-propyl-pyridin-3-yl]ethanol 




FP005ME 
FP007SE 



Merck 
Serono 



1-(2-Trlfluoromethoxyptienyl)-2-nitroetiianone 




FP013NO Novartis 



FP014SC Sctiering Tetraettiyl[(3-tiydroxy-2-pyridyl)-amino]-mettianediptiosptionate 





100 20 Liver, kidney 4-Phenylpyridine Glucagon receptor antagonist 



350 15 Liver Carbonylnitrometfiane Carbonylnitromettianes 

1000 100 Kidney liver 



Cyclosporin derivative - 



1120 280 Liver Bisphosphonate Inhibitor of osteoclast activity 
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Table S2 Primer sequences 

Sequences of primers used to generate Dio3 and Actb probes used in qRT-PCR and PGR 
experiments. 

Gene Experiment Direction Sequence (5'^3') 



Dio3 


qRT-PCR 


Forward 


AGAGTGGCACCATCATGTACCA 






Reverse 


CCAAGTGCGCAACTCAGACA 


Dio3 


PGR and sequencing 


Forward 


GTGGTCGGAGAAGGTGAA 






Reverse 


AGAGGCGCTCAAAATAGG 


Actb 


qRT-PCR 


Forward 


CGTGAAAAGATGACGCAGATCA 






Reverse 


CACAGCCTGGATGGCTACGT 


Actb 


PGR 


Forward 


ATGGTGGGTATGGGTCAGAA 






Reverse 


AGATGTGCTGGAAGGTGGAC 
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Table S3 mRNAs changed in translational activity following treatment with PredTox compound FP014SC 

Differentially translated mRNAs were identified following sucrose density fractionation and polysome profiling of liver samples treated witfi FP014SC. Those mRNAs ttiat showed a statistically 
significant (P < 0.05) difference in ribosomal occupancy in at least one subpool of fractions when compared with vehicle-treated samples are listed. Significant values are denoted by bold type. 
Dio3 is highlighted. Fold changes are colour-coded as indicated in the key 
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Decreased translation of Dio3 mRNA following liver damage 



Table S4 Canonical pathways that were significantly perturbed (two or more genes) following treatment with PredTox compound FP014SC 

All mRNAs that showed a change In ribosomal occupancy across the gradient following polysome profiling were uploaded for pathway analysis to Ingenuity® IPA software. Those pathways under 
the control of RXR are In bold (FXR, farnesold X receptor; IL-1 , Interleukin 1 ; LPS, llpopolysaccharlde; LXR, liver X receptor; PPARa, peroxisome-prollferator-actlvated receptor a; PXR, pregnane X 
receptor). 



Ingenuity® canonical pathway 


- log (P value) 


Nui 


Fatty acid metabolism 


8.71 


15 


LPS/IL-1 mediated inhibition of RXR function 


2.67 


10 


Production of nitric oxide and reactive oxygen species in macrophages 


3.10 


10 


Clathrin-mediated endocytosis signalling 


2.99 


10 


G-protein-coupled receptor signalling 


0.319 


9 


mTOR signalling 


2.47 


9 


Axonal guidance signalling 


0.689 


9 


Glucocorticoid receptor signalling 


1.51 


9 


Tryptophan metabolism 


2.81 


8 


Protein Idnase A signalling 


0.834 


8 


PPARa/RXRa activation 


2.11 


8 


Metabolism ot xenobiotics by cytochrome P450 


3.18 


8 


LXR/RXR activation 


3.15 


8 


RAR activation 


2.11 


8 


elF2 signalling 


1.94 


8 


Regulation of elF4 and p70™ signalling 


2.57 


8 


Protein ubiquitination pathway 


1.16 


8 


PI3K signalling in B-lymphocytes 


2.14 


7 


Role of macrophages, fibroblasts and endothelial cells in rheumatoid arthritis 


0.613 


7 


Molecular mechanisms of cancer 


0.475 


7 


Calcium signalling 


1.58 


7 


Xenobiotic metabolism signalling 


0.718 


7 


Dopamine-DARPP32 feedback in cAMP signalling 


1.71 


7 


ERK/MAPK signalling 


1.47 


7 


Acute-phase response signalling 


1.67 


7 


Valine, leucine and isoleucine degradation 


3.11 


6 


GREB signalling in neurons 


1.17 


6 


B-cell receptor signalling 


1.48 


6 


Breast cancer regulation by Stathmini 


0.981 


6 


Wnt ;6-catenin signalling 


1.19 


6 


Bile acid biosynthesis 


3.58 


6 


Arachidonic acid metabolism 


1.75 


6 


RhoGDI signalling 


0.773 


5 


Ephrin receptor signalling 


0.794 


5 


Huntington's disease signalling 


0.530 


5 


Purine metabolism 


0.330 


5 


Linoleicacid metabolism 


2.07 


5 


Glycolysis/gluconeogenesis 


1.81 


5 


Fatty acid elongation in mitochondria 


4.93 


5 


Phospholipase G signalling 


0.468 


5 


Atherosclerosis signalling 


1.33 


5 


IL-12 signalling and production in macrophages 


1.15 


5 


Actio cytoskeleton signalling 


0.554 


5 


Sertoli cell-Sertoli cell junction signalling 


0.745 


5 


Semaphorin signalling in neurons 


1.97 


4 


NRF2-mediated oxidative stress response 


0.457 


4 


P2Y purigenic receptor signalling pathway 


0.848 


4 


Tyrosine metabolism 


1.46 


4 


Pyruvate metabolism 


1.62 


4 


Human embryonic stem cell plurlpotency 


0.753 


4 


GNRH signalling 


0.821 


4 


Role of NFAT In cardiac hypertrophy 


0.439 


4 


Synaptic long-term potentiation 


1.02 


4 


Colorectal cancer metastasis signalling 


0.250 


4 


Glioblastoma multiforme signalling 


0.623 


4 


PI3K/Akt signalling 


0.867 


4 


PTEN signalling 


0.979 


4 


Hereditary breast cancer signalling 


0.886 


4 


Tight junction signalling 


0.585 


4 


IL-8 signalling 


0.491 


4 


Glycerolipid metabolism 


1.05 


4 


FXR/RXR activation 


1.25 


4 


Inositol phosphate metabolism 


0.616 


4 


AMPK signalling 


0.737 


4 


Mitochondrial dysfunction 


0.737 


4 


Type 2 diabetes mellitus signalling 


0.906 


4 



of altered genes 
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K. M. Dudek and others 



Table S4 Continued 



Ingenuity® canonical pathway - log (P value) Number of altered genes 



Ceramide signalling 1.39 4 

Aldosterone signalling in epithelial cells 0.579 4 

Amyotrophic lateral sclerosis signalling 1.08 4 

Signalling by Rho family GTPases 0.266 4 

Sphingosinel-phosphate signalling 0.926 4 

Germ cell-Serfoli cell junction signalling 0.573 4 

p38 MARK signalling 0.590 3 

Circadian rhythm signalling 1.84 3 

FGF signalling 0.787 3 

Pentose and glucuronate interconversions 1.55 3 

GM-CSF signalling 1.08 3 

Neuropathic pain signalling in dorsal horn neurons 0.605 3 

iCOS-iCOSL signalling in T-helper cells 0.590 3 

PKCe signalling in T-iymphocytes 0.526 3 

Prostate cancer signalling 0.820 3 

Cyclins and cell-cycle regulation 0.843 3 

Aryl hydrocarbon receptor signalling 0.401 3 

Small-cell lung cancer signalling 0.944 3 

Thyroid cancer signalling 1.44 3 

Non-small-cell lung cancer signalling 1.03 3 

Retinol metabolism 1.44 3 

TR/RXR activation 0.777 3 

Galactose metabolism 1.47 3 

Alanine and aspartate metabolism 1.55 3 

Dendritic cell maturation 0.304 3 

0X40 signalling pathway 0.766 3 

Cdc42 signalling 0.293 3 

Type 1 diabetes mellitus signalling 0.567 3 

Insulin receptor signalling 0.437 3 

Natural killer cell signalling 0.605 3 

Oxidative phosphorylation 0.368 3 

Role of 0ct4 in mammalian embryonic stem cell pluripotency 1.42 3 

Cardiac ,6-adrenergic signalling 0.426 3 

fMLP signalling in neutrophils 0.575 3 

OCRS signalling in eosinophils 0.533 3 

CXCR4 signalling 0.330 3 

Paxillin signalling 0.629 3 

LPS-sfimulafed MARK signalling 0.917 3 

Regulation of actin-based motility by Rho 0.787 3 

Rac signalling 0.605 3 

p70^ signalling 0.476 3 

NGF signalling 0.605 3 

Androgen and oestrogen metabolism 0.832 3 

Sulfur metabolism 2.38 3 

Fructose and mannose metabolism 1.30 3 

Renal cell carcinoma signalling 0.958 3 

Hypoxia signalling in the cardiovascular system 1.00 3 

RhoA signalling 0.258 2 

Neurotrophin/Trk signalling 0.498 2 

ATM signalling 0.671 2 

Arginine and proline metabolism 0.466 2 

Ubiquinone biosynthesis 0.490 2 

Basal cell carcinoma signalling 0.490 2 

Factors promoting cardiogenesis in vertebrates 0.372 2 

Role of Nanog in mammalian embryonic stem cell pluripotency 0.266 2 

Glioma signalling 0.329 2 

Melatonin signalling 0.474 2 

Urea cycle and metabolism of amino groups 1.01 2 

Oestrogen-dependent breast cancer signalling 0.581 2 

p53 signalling 0.339 2 

Bladder cancer signalling 0.361 2 

Androgen signalling 0.254 2 

Cell cycle: Gi/S checkpoint regulation 0.624 2 

Chronic myeloid leukaemia signalling 0.334 2 

Cell cycle regulation by BTG family proteins 0.942 2 

HGF signalling 0.309 2 

Pancreatic adenocarcinoma signalling 0.278 2 

PXR/RXR activation 0.498 2 

IL-1 signalling 0.344 2 

Glutamate receptor signalling 0.613 2 
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Decreased translation of Dio3 mRNA following liver damage 



Table S4 Continued 



Ingenuity® canonical pathway - log (P value) Number of altered genes 



Butanoate metabolisrr 


0.581 


2 


Telomerase signalling 


0.305 


2 


Cross-talk between dendritic cells and natural killer cells 


0.355 


2 


Cytotoxic T-iymphocyte-medlated apoptosis of target cells 


0.409 


2 


Allograft rejection signalling 


0.397 


2 


Caveoiae-mediated endocytosis signalling 


0.458 


2 


Virus entry via endocytic patliways 


0.355 


2 


Neuroprotective rale of TH0P1 in Alzheimer's disease 


0.848 


2 


FcsRi signalling 


0.291 


2 


PDGF signalling 


0.482 


2 


CD40 signalling 


0.552 


2 


TNFR1 signalling 


0.724 


2 


Angiopoietin signalling 


0.533 


2 


Parkinson's disease signalling 


1.45 


2 


Apoptosis signalling 


0.361 


2 


Retinoic acid-mediated apoptosis signalling 


0.552 


2 


Citrate cycle 


1.08 


2 


Activation of IRF by cytosollc pattern recognition receptors 


0.552 


2 


Dopamine receptor signalling 


0.444 


2 


CDK5 signalling 


0.372 


2 


Regulation of iL-2 expression In activated and anergic T-lymphocytes 


0.423 


2 


CD28 signalling in T-heiper ceils 


0.240 


2 


Agrin Interactions at neuromuscular junctions 


0.507 


2 


FAK signalling 


0.372 


2 


PAK signalling 


0.361 


2 


Renin-angiotensin signalling 


0.266 


2 


Role of BRCA1 in DNA damage response 


0.624 


2 


Sphingolipid metabolism 


0.397 


2 


Keratan sulfate biosynttiesis 


0.635 


2 


Inositol metabolism 


2.41 


2 


G-protein signalling mediated by Tubby 


0.963 


2 


14-3-3-mediated signalling 


0.226 


2 


Starch and sucrose metabolism 


0.507 


2 
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